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Abstract 

The classical theory of radiating point-charges is revisited: the retarded 
potentials, fields, and currents are defined as nonlinear generalized functions 
and all calculations are made in a Colombeau algebra. The total rate of 
energy-momentum radiated by an arbitrarily moving relativistic point-charge 
under the effect of its own field is shown to be rigorously equal to minus 
the self-interaction force due to that field. This solves, without changing 
anything in Maxwell's theory, numerous long-standing problems going back 
to more than a century. As an immediate application an unambiguous 
derivation of the Lorentz-Dirac equation of motion is given, and the origin of 
the problem with the Schott term is explained: it was due to the fact that the 
correct self-energy of a point charge is not the Coulomb self-energy, but an 
integral over a delta-squared function which yields a finite contribution to the 
Schott term that is either absent or incorrect in the customary formulations. 

03.50.De Classical electromagnetism, Maxwell equations 

1 Introduction 



Gauss's theorem in its various forms is a major tool in field theory because it 
enables to transform integrals over closed surfaces into integrals over the volume 
contained within them. This enables, for example, to express Maxwell's equations 
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in integral rather than differential form, and thus to establish precise relations 
between the total flux of fundamental field quantities through closed surfaces and 
the volume integrals of quantities that are differentially related to them. 

As long as all field quantities are smooth enough, and provided all surfaces are 
sufficiently regular and simple, there are no problems: Gauss's theorem and its 
numerous corollaries can be applied without restriction and with no great technical 
difficulties. But this is not the case when there are singularities in the fields. 

Consider for example the flux of energy-momentum radiated by a point-charge 
under the effect of its own field through a 3-surface enclosing its world-line, and 
the related 4- volume integral of the self-interaction force on this point-charge due 
to that field: according to Gauss's theorem they should be equal — provided all 
problems related to singularities are properly taken care of since the fields are 
infinite on the world-line. Indeed, let T^ u by the symmetric energy-momentum 
tensor, and 3 S the 3-surface enclosing a 4-volume 4 f2 surrounding the world-line 
between the proper-times n and r 2 . The radiated energy-momentum flux is then 

P M (£):=//7 TV = ^ dr [[[ d 3 Q (1.1) 

lA/3 S J T1 JJJ3Q 



where Gauss's theorem was used to write the right-hand side. Using the definition 

ofT M „ 

= -^-{F^F» - -gTF^F^), (1.2) 

and Maxwell's equations 



^ + ^ + ^f=0, and d»F, u = -4n.J„ (1.3) 

the right-hand side of Eq. (|l.ll) becomes 

/ 2 dr fff d 3 n (-F^r) := f 2 dr P^Q) = F M (fl). (1.4) 
Jti JJJ^n Jt x 

Therefore, comparing with the left-hand side of Eq. (11.11) . Gauss's theorem implies 

P M (£) = P^) = -Q M (Q), (1.5) 

where 

Q„(Q) := JJJ (FtoF^r. (1.6) 

Equation Eq. (11.51) states that the rate of energy-momentum flux through the 
3-surface 3 E at the time r, calculated as the proper-time derivatives of either the 
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3-surface integral -P M (S) or of the 4-volume integral P M (f2), is equal to minus the 
integral Q^Q) of the self-interaction force density F IMU J U over the 3-volume 3 fi 
enclosed by the 3-surface at that time. 

However, when P M (£) and -P^(^) are evaluated using for F M „ the customary 
Lienard-Wiechert field strengths and for J M the customary current density of an 
arbitrarily moving relativistic point-charge, what is not without difficulties and 
ambiguities since both F^ v and J u are unbounded on the world-line, one finds that 
they are not equal. 

For example, following Lorentz [1J, the self-force is usually calculated non- 
relativistically and then generalized to the relativistic case. This gives 

P„(n) = lim ~Z^ - \e 2 A 2 Z^ - |e%, (1.7) 

where the term containing the square of the acceleration A 2 = Z^Z^, known as the 
'Larmor term,' is the invariant radiation rate, and the term in the third proper-time 
derivative of the position, Lorentz's frictional force, is known as the 'Schott term.' 

On the other hand, the energy-momentum flux rate can be found by either 
generalizing a non-relativistic calculation, or, following Dirac ED, by a direct 
covariant 4-dimensional calculation. This gives 

P M (E) = gm~£„ - 2 -e 2 A 2 Z^ (1.8) 

in which the Schott term is missing and where the diverging quantity e 2 /2£ is 
called the 'Coulomb self-energy' or 'electrostatic self-mass' of the point-charge, 
which differs from the 'electromagnetic self-mass' 2e 2 /3£ which appears in (11.71) 
by a factor 4/3. 

Apart from this notorious '4/3 problem,' the main difference between P^E) 
and P M (f2) is that Pj^O) is a Minkowski force, that is a genuine relativistic force 
satisfying the orthogonality condition 

P^ M = 0, (1.9) 

on account of the kinematical identities Z^Z^ = 1, Z^Z^ = 0, and z\ t Z^ = —A 2 , 
whereas P M (S) is not such a force because there is no Schott term. In fact, this 
missing term is a major problem because, as noted by Dirac, the derivation of 
(|1.8I) is quite simple, and there is no way to get the Schott term in this calculation 
without making a supplementary assumption. 

These problems, as well as other inconsistencies, have led to hundreds of 
papers and endless controversies which we are not going to review and discuss 
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here. Indeed, while theses difficulties are known, there is a large consensus that 
the finite part of the force (11.71) is the correct relativistic form of the self-force 
on an arbitrarily moving point-charge, i.e., the so-called 'Abraham' or 'von Laue' 
self-interaction force |e 2 (Z M + A 2 Z^). Similarly, there is a large consensus that 
the divergent terms in (11.81) and (11.71) should be discarded or 'renormalized', and 
that the equation motion of a point-charge in an external electromagnetic field F™ 4 
is 

mZ, = \e 2 (z\ + A%) + eF^Z\ (1.10) 

i.e., the Abraham-Lorentz-Dirac equation of motion. 

In this paper we show that the fundamental reason for discrepancies such as 
between (11.81) and (11.71) is mathematical, and due to two shortcomings of the 
customary calculations: 

• incomplete characterization of the singularities of the fields and current at 
the position of the point charge — everything should be defined not just as 
distributions but as nonlinear generalized functions EH); 

• use of formalisms such as elementary calculus or Schwartz distribution 
theory in which the product of distributions is not defined — this is a major 
defect since the energy-momentum tensor (|1.2I) is nonlinear in the fields so 
that all calculations should be made in a Colombeau algebra [EL IS 131 [HI - 

Once all singularities are properly defined, what basically consists of redefining 
the Lienard-Wiechert potential as a Colombeau nonlinear generalized function 
BH QZ5), the calculations of the self-energy-momentum radiation rate and of the 
self-interaction force are straightforward (but technically complicated) and the 
equality (11.51) is obtained for a strictly point-like electric charge in fully arbitrary 
relativistic motion without making any modification to Maxwell's theory. 

The importance of this result lies is the fact that it leads to the solution of 
a number of difficulties which have often been interpreted as suggesting that 
Maxwell's electrodynamics of point-like electrons is not fully consistent, and 
required to be modified at the fundamental level. To illustrate that this is not 
necessary, it will be shown how the Lorentz-Dirac equation of motion derives 
from our rigorously calculated self-interaction force. 

Because the full calculation leading to the proof of the equality (11.51) for a 
relativistic point-charge is very long and complicated, the details are given in 
another paper IfTTTl . Nevertheless the general physical concepts, the principle of 
the mathematical methods used, and the main steps in the calculation are explained 
in the present paper. 
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2 Definitions of the fields and of the geometry 



Let Z^ be the 4-position of an arbitrarily moving relativistic point-charge, and 
a point of observation. At the point X M the 4-potential A^, field F^ v , and 4-current 
density J M are functions of the null interval i? M between X M and Y^, i.e., 

:= - Z M , such that R^R? = 0, (2.1) 

as well as of the 4- velocity Z^, 4-acceleration Z^, and 4-biacceleration Z^ of the 
charge, to which three invariants are associated: Z^R^, Z^R^, and Z^R^. The 
first one is called the retarded distance, 

£ := Z,R», (2.2) 

which enables to introduce a 'unit' null 4-vector defined as 

K^(6,<f>) :=V6 (2-3) 
and the so-called acceleration and biacceleration invariants defined as 

and x := Z^K". (2.4) 

The three spatial variables {£, 0, 0} and the proper-time r define a causal 
coordinate system that is most appropriate for expressing the fields of a point- 
charge moving according to the world-line prescribed by Z^(t). As explained in 
Refs. [9] and IfTOl . the retarded potential of such a point-charge must be written as 
a nonlinear generalized function, i.e., 



4. = e y T (0 



(2.5) 



where T(£), which is absent in the customary Lienard-Wiechert formulation, is 
the generalized function defined as 

f undefined r < 0, 
r = 0, and 4 T ( r ) = <K r )> (2-6) 

ar 

r > 0, 



which explicitly specifies how to consistently differentiate at £ = 0. Instead of the 
customary Lienard-Wiechert field, the field strength is then |[9l 



F,u = e(^T(0 + (1 - <)^(T(£) - £5(0) 
A ' " Z "-T(0 - (1 - kO^(T(0 - £5(0) 



(2.7) 
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which apart from the presence of the T-function factors, has additional 5-like 
contributions. As a consequence of these T and 5 factors, the current density 
deriving from this field, i.e., 



is significantly more complicated than the customary current density, but turns 
out to be locally conserved, whereas the customary one is not [|9l [T0]|. In this 
equation, as in Eqs. (|2.6I) and (12.71) . the condition r = r r implies that all quantities 
are evaluated at the retarded proper time r r . In the following, for simplicity, this 
condition will be specified explicitly only for the main equations. 

Having defined the field F^ v and the current density J M to be used to calculate 
the integrands of Eq. (11.11) . (11.41) and (11.61) . it remains to define the geometry of a 
3-surface enclosing the world-line in order to calculates the integrals themselves. 

In view of this we begin by ignoring the fact that is singular on the world- 
line. We are then free to chose any closed 3-surface intersecting the world-line at 
two fixed proper times n and r 2 . This is because by Gauss's theorem the field- 
integrals over the 4-volume enclosed between any two such 3-surfaces will give 
zero since d v F {W = 0, and thus d v T iiV = 0, in that 4-volume. To be consistent 
with the already defined causal coordinate system {r, £, 8, (f\ we therefore take a 
3-surface consisting of one 3-tube of constant retarded distance £ = Cst between 
the times n and r 2 , and, to close that 3-tube at these times, two 3-spheres r = Cst 
at the times n and r 2 between the world-line and the radius of the 3-tube. We 
therefore do not take a light-cone (defined as r r := t — £ = Cst) to close the tube 
at both ends as is often done to simplify calculations when only surface integrals 
are calculated: such a choice would make the corresponding 4-volume integrals 
significantly more complicated. Besides, as the formulation of our fields involves 
T and 5 functions whose argument is £, taking surfaces and volumes different 
from those defined by constant values of the causal coordinates £ and r would 
make the evaluation of the fields singularities at £ = even more complicated 
than they already turn out to be. Finally, with our choice, the 4-volume element is 
the simplest possible, namely 



so that the 4-volume integrals are simply proper-time integrals times ordinary 

3- volume integrals with the radial variable equal to £. 

Now, to deal with the singularities on the world-line, we replace the just defined 

4- volume, that is such that £ varies between £ = and some maximum £ = £ 2 , 
by a 'hollow' 4-volume in which £ varies between £ = £i 7^ and the maximum 




(2.8) 



= dr d<p d0sin(0) £ 2 , 



(2.9) 
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£ = £ 2 - This 4- volume will then be comprised between two 3-tubes of radii ^ and 
£2, where £1 will eventually tend to at the end of the calculation. Thus, we will 
calculate the 4- volume integral 

dr /// d 3 tt(...) = dr d^ 2 // du (...), (2.10) 

T\ JJJ J Tl «/ ^1 */«/ 

in which (for practical reasons) the integration over the ordinary solid angle duo = 
d<p dO sin(6 l ) will be done first in general. 



In this geometry, the normal to the 3-tube £ = Cst is 

Nl = (1 - kOK, + Z^ 



5 

r=r r 



(2.11) 



and the integrals over the internal and external 3-tubes are 

d^l (...) = dr J! du;£ 2 iVj(-)|^ 



(2.12) 



where £1 is the internal radius of the tube that is sent to zero at the end of the 
calculation. 



Similarly, the normal to the 3-sphere r = Cst is 



5 

T = T r 



(2.13) 



and the combined integrals over the two hollow 3-spheres enclosing the hollow 
4-volume at both ends can be written 



(2.14) 



3 Mathematical methodology 



The proof of Eq. (11.51) involves calculating products of the generalized functions 
T, S, and 5' such as T 2 , ST, S'T, 5'5, and S 2 , possibly multiplied by £ n with 
n G Z. If these generalized functions were defined as Schwartz distributions these 
products would not be defined and all sorts of inconsistencies would arise. But as 
we define them as nonlinear generalized functions we can use Colombeau's theory 
which enables to multiply, differentiate, and integrate them as if we were working 
with ordinary smooth functions 
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In Colombeau's theory the generalized functions T, 5, etc., are embedded in 
an algebra by convolution with a 'mollifier function' r)(x), i.e., an unspecified real 
smooth function whose minimal properties are given by 

J dz r](z) = 1, and J dz z n r}{z) = 0, Vn = 1, q E N. (3.1) 

This function can be seen as describing details of the singularities at the infinites- 
imal level. These details are possibly unknown, but not essential for getting 
meaningful results which are generally obtained by integrations in which only 
global properties such as (|3.1I) are important. However, as will be the case in 
this paper, further constraints can be put on the mollifier, for instance to meet 
requirements imposed by the physical problem investigated. 

Of course, working in a Colombeau algebra requires some precautions: all 
'infinitesimal' terms aring as a result of a differentiation or any other operation must 
be kept until the end of the calculation — contrary to calculating with distributions 
where such terms may possibly be discarded. Similarly, distributional identities 
like S'(x) = —6(x)/x, which was used to write the conserved Lienard-Wiechert 
current density in the form (|2.8I) . have to be avoided to remain as general as 
possible. We will therefore use J M in the form ifTOl Eq. (4.17)] 

= ^(-(! " 2"Z)zA6' + {Z„ - ixK^zd + kN^(U - 5')), (3.2) 

where iVj is the normal to the 3-tube £ = Cst, i.e., Eq. (12.111) . 

To calculate the integrals (12.101) . (12.121) . and (12.141) . we introduce a notation 
that will be systematically used in the following. That is, we rewrite the field 
strength as a sum of two terms 

F, u (r,0 = F^(r,Or(0 + K^Om, (3-3) 

where the upper indices ' T ' and 1<5 ' refer to the T part (i.e., the customary Lienard- 
Wiechert field) and respectively to the 5 part of Eq. (12.71) . The energy-momentum 
tensor T^ v consists then of four terms, and consequently the surface integrals 
(12121) and (12T41) comprise also four terms, i.e., Pj ,T (S), Pj ,<5 (£), ^' T (£), and 
P / f ,<5 (S), in which the factors T 2 , T5, ST, and 5 2 appear. Similarly, the volume 
integral (|2.10l) involves then two terms Pj' s (£l) and P^' 5 (Cl), in which the factors 
ST and S'T, and respectively 55' and 5 2 , appear. 

As for the integrations themselves, we perform them according to a similar 
pattern, in which all £ integrations are left for the end because both 5 and T are 
functions of £ which are singular as £ — » 0. Thus we begin by making the angular 
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integrations, which together with the complicated algebra implied by the products 
of the many field components is the most laborious part of the calculation. This 
work is however greatly facilitated by the biquaternion techniques due to Paul 
Weiss explained in Refs. lITOl and |fTTfl. which take maximum advantage of the 
cancellations between null 4-vectors, and which allow a maximal separation of 
variables and factorization of all quantities in terms of factors depending solely of 
kinematical or angular variables [12J. The resulting surface and volume integrals 
are given in Secs.|4]and[H 

Next we transform the resulting surface integrals in such a way that they have a 
form directly comparable to the volume integrals, namely a double integration over 
£ and r. This enables, after dropping the r integration, to compare the remaining 
^-integrals, which should be equal according to Eq. (11.51) . The integrands of 
Eqs. (12.121) and (12.141) . which are initially integrated over r and respectively £ 
alone, have therefore to be transformed and combined in non-trivial ways. This is 
made possible by the fact that all kinematical quantities are implicit functions of 
the retarded time r r = r — £, so that their r and £ derivatives are related as 

— E{r r ) = -Q-Mr r ), or E{r r ) = -E{r r )', (3.4) 

where the dot and the prime denote differentiation with respect to r and to £, 
respectively, and where E(r r ) is any expression function of the variable r r . A 
particularly interesting intermediate result in the transformation of the surface 
integrals is given as Eqs. (14.7144.91) . 

Finally, the volume integrals given in Sec. [5] are transformed in such a way 
that there is no 5' anymore. This is done using the identities (5 2 )' = 255' and 
(f 5T)' = 5T + £5'T + £5 2 , which enable to integrate (E2) and (O by parts to 
eliminate 5'. 

At this stage all expressions are comparable. It is then found that Eq. (11.11) 
holds exactly as it should according to Gauss's theorem, and that both -P M (S) and 
P M (fi) can be written in the form (I6.3H6.6I) . Thus Eq. (11.51) holds as well and 
P^E) = can be written in the form (I7.H - I7.4I) . However, since we have 

not yet made the final £ integrations, these equalities are only a check that all 
algebraic calculations and angular integrations are correct, and that the integrands 
have been properly transformed using Eq. (13.41) and the rules of partial integration. 
The precautions taken to properly deal with the singularities, and the postulate 
that everything is calculated in a Colombeau algebra, will become essential when 
making these £ integrations and taking the limits £ 2 — ► oo and £x — > 0, what will 
be done in Sec. [7] 

Finally, until Sec.|71 we will take e = 1 for the electric charge to avoid the 
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factor e 2 which is multiplying all surface and volume integrals. 



4 Surface integral: the energy-momentum flow 



The four contributions Pj' r , Pj' 5 , P^,' r , an d P^ 5 to the surface integral -P M (S), 
i.e., Eq. (11.11) or the proper-time integral of the left-hand side of (11.51) . are: 



T2 dr (i^-l^ lT 3 



+ 



71 



1 



(4.1) 
(4.2) 



+ / ^ 



T2 



(4.3) 
(4.4) 



x2£ 2 ~ At 6£~" 6'" M 10 



(4.5) 
(4.6) 



When £ — > the expression inside the parentheses in (14.11) is equal to (11.81) . i.e., 
the usual energy-momentum flow rate through a 3-tube surrounding the world-line 
El [12l . This illustrates how much more complicated the rigorous formalism and 
calculations are than the customary ones. 

These integrals can be transformed according to the methods outlined in Sec. [3] 
and, in particular, (|4. 1114.41) can be written 



P, T ' T fS) 



T-2 



dr 



d^[+^-^A 2 Z u )26T, 



and 



^ t ' 5 (s) + p m 5 ' t (s; 



— Za — —Zn )2(5T 

2£ M 6 M/ 



(4.7) 

(4.8) 
(4.9) 
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When adding (14.71) and (14.81) a remarkable cancellation takes place: the l/2£ 
diverging terms cancel each other. This means that the Coulomb self-energy does 
not contribute to the total surface integral Indeed, the diverging term which 

will be interpreted as the 'self-mass' will not be the usual Coulomb self-energy, 
i.e., the first term in (14.11) . but the leading divergence in (14.61) . which is an integral 
over a <5 2 -function. This is consistent with what was found in Ref. Q for a point- 
charge at rest, namely that its self-energy is not the usual Coulomb self-energy 
but an integral over a 5 2 -function, so that the self-mass is fully concentrated at the 
position of the singularity. 



5 Volume integral: minus the r-integrated self-force 

The two contributions Pj' s , and P^' 5 to the volume integral (11.41) . i.e., minus the 
proper-time integral of the right-hand side of (11.51) . are: 

Pj> s (Q) = £ dr£ di [—Z, - ^A%)ST (5.1) 

dr J di [--A 2 Z^8'T, (5.2) 

and 

P 5 /m = jT dr J' d£ + ^A 2 Z, + ^A 2 Z^e55' (5.3) 

+ f* dr J 2 di + \A 2 Z, + ^A 2 Z, - £(ZuZ")Zp)t#. 

(5.4) 

These integrals are rather complicated, what is of course necessary to match the 
complexity of the surface integrals. In fact, this could be expected from the 
conclusion of Ref. 13, where is was explained that the self-force density on an 
arbitrarily moving point-charge is given by F iW 3 v where is the complete 
retarded field (|2.7I) . and where J u is the full locally-conserved current density 
(EE) or (Q. 
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6 The general energy-momentum flow 



By fully transforming the surface and volume integrals according to the methods 
sketched in Sec.[3l the sums 

P„(E) = P M T < T (£) + P^(E) + P*> T (£) + If (E), (6.1) 
P,(Q)=P^ 5 (Q) + P S /(Q), (6.2) 

can be shown to by equal, i.e., P M (E) = P M (0) := P M as required by (11.11) . and 
can both be put in the form: 

P,(n, r 2 , 6,6) = jH dr J 2 d£ (-^'m - - |(^ 2 ^)')5T (6.3) 



n 

T2 



1 ~ 1 „2~ £ 



6 



ft 



1 ... £ 2 • 



(6.4) 
(6.5) 



/ T2 dr I 2 di [Lz - \z» + - £(A%y)&, 



n 



(6.6) 

where all kinematical quantities are evaluated at the retarded proper time r r = r — £. 

For convenience we refer to (I6.3H6.6I) as the 'general energy-momentum flow' 
even though we could change its overall sign and call it the 'general r-integrated 
self-force.' By 'general,' we mean (i) that this flow is a function of the parameters 
r i, r 2,£,i, and £ 2 ; an d (ii) that the function T, from which 5 derives by differentia- 
tion, belongs to a class of the nonlinear generalized function which may be further 
specified by physical conditions leading, for instance, to additional constraints on 
the Colombeau mollifier r\. 



7 The general self-interaction force 

By differentiating P^(ti, t 2 , £i, £2) with respect to the proper time, that is by 
suppressing the r integrations in (16.3H6.6I) . we get the equality (11.51) . i.e., P M (E) = 



12 



P^Q) = -Q„ with 



1 - 



e 



10" 



+ e z 



5T 



di e 2 



Si 



2 M 



6 



e 



6 (A 2 z„y + ^(X'z»y)? 



(7.1) 
(7.2) 
(7.3) 

(7.4) 



where all kinematical quantities are evaluated at the retarded proper time r r = r — £, 
and where in view of applications we have we have made the e 2 factor explicit. 

Equation (17. 1117.41) is the main result of this paper: it is an exact closed-form 
expression giving, in the limits £ 2 — > oo and £ x — > 0, the total self-interaction 
force on an arbitrarily moving relativistic point-charge. Because it leads to the 
Abraham- von Laue self-interaction force when taking these limits we call it the 
'general self-interaction force.' 

To calculate the limits £ 2 — > oo and £j — > in (|7. 1147 .41) we first remark 
that if Z(r) is indefinitely differentiable in r, which implies that the velocity, the 
acceleration, and the biacceleration are smooth functions of r, the kinematical 
expressions inside the big parentheses are all smooth as functions of r r , that is as 
functions of £ as well. Then, because 



£5T 



0. 



and 



0. 



(7.5) 



the boundary terms (|7.11 - I7.2I) give zero. As for (17 .3147 .41) we use the formulas dH|7]| 



dr T(r)5(r)T(r) = -T(0), 



dr S 2 (r)T(r) 



iimq ]^ + qi]T'(o), 



where Cr i and Cm are functions of the Colombeau mollifier r], i.e., 

" + 0O 



Cu 



dx 7]" 



-x) 



and 



[i] 



dx xi] (—x). 



(7.6) 
(7.7) 

(7.8) 



In (17 .6117 .71) the function T(r) is a smooth function with compact support, which 
implies that T(±oo) = 0, what is the case of the kinematical expressions in (17.31 - 
17.41) if we suppose that Z(±oo) = 0, i.e., that the velocity is zero at r = ±oo. 
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Applying (17.54 - 47.71) to (I7.H - I7.4I) we finally get for the total self-interaction 
force Q^ elf := Q^O, +00) the expression 

Qf (r) = e 2 (- lim ^4( T ) + %^( T ) + l^(r) + ^ 2 (r)Z,(r)) , (7.9) 

where all kinematical quantities are evaluated at £ = 0, i.e., on the world-line 
at the proper-time r, and where we used (13.41) to replace the ^-differentiation in 
(T7.7I) by a r-differentiation to get the term +|C[i]Z /U (t). We therefore see that 
the <5 2 -integral (17.41) contributes two terms to the self-interaction-force (T7.9I) : a 
diverging contribution which can be interpreted as a self-mass as in the customary 
calculations, and a finite additional contribution to the Schott term which does 
not arise in the customary calculations. If these contributions do not need to be 
explicitly shown, one can rewrite (17.91 ) using (17.71) in the more concise form 

Qf(r) = e 2 (-- J 8%)Z,{r r ) df+ -Z„(r) + ~-4 2 (t)^(t)) , (7.10) 
which like (I6.3H6.6I) is valid in any algebra of nonlinear generalized functions. 



8 The Lorentz-Dirac equation of motion 

An immediate application of the general formalism developed in this paper and 
of the self-interaction force <5^ elf is the derivation of the equation of motion of a 
point-particle in arbitrary motion in an external electromagnetic field F^. Indeed, 
the total force on a point-particle under the effect of its own field (12.71) and of an 
external field is 

m Z,{r) = JJJ d 3 n (F, u + F^)r = Qf(r) + eF^(Z)Z^r), (8.1) 
with the subsidiary condition 

Z^ = 1, (8.2) 

which specifies that the world-line Z(r) is the world-line of a point in Minkowski 
space. In writing (18.11) we postulated that the total force could be written m Z^ 
where m is a constant called the 'non-renormalized' mass, and that the external 
field Fj^(X) is a smooth function of X — or at least smooth at any point X 
near or on the world-line Z{f) — so that the volume integral of the external force 
density, i.e., F^(X)J U with J u given by (HH) or (O, reduced to eF^{Z)Z v 
after an elementary integration over a three-dimensional 5-function. 
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Substituting Q self from (17.91 ) we get therefore 

(m + lim^)^ = e 2 (^Z, + + ^ 2 2,) + eF-^, (8.3) 

where we have put the diverging self-mass term on the left-hand side. We now 
use the subsidiary condition condition (18.21) . which by differentiation leads to the 
identities Z^ = and Z^ = -A 2 . Thus, contracting dOl) with Z^ and 
using the antisymmetry of the the external field we obtain 

C m 1 2 

-£!■ + - = -. (8.4) 
2 6 3 

This equation is satisfied if the Colombeau mollifier 77 is such that the moment 

P+OO 

Cm = / dx xrj 2 (— x) = 1, (8.5) 



which is an additional condition supplementing the minimal conditions (13.11) . 
Consequently, Eq. 08.31) becomes 

(mo + lim ^Z„ = 2 -e\z\ + A 2 Z^) + eF^Z\ (8.6) 

which after substituting the 'renormalized' or 'experimental' mass 

C 

m := m + lim — — , (8.7) 

e-^o 2e 

is known as the Lorentz-Dirac equation of motion, i.e., Eq. (11.101) . where everything 
is evaluated at the proper time r. 

To fully appreciate what has been achieved, let us compare our derivation to a 
particularly elegant approach due to Barut |[T3l . The idea is to calculate the self- 
interaction force by analytically continuing the Lorentz force of a point-charge on 
itself, i.e., eF^ f Z u (r) where Fpf is the customary Lienard-Wiechert field, to a 
point X = Z(t + e) close to the world-line. If this done, one easily finds lfT3l 
Eq.8] 



QSLCO = e 2 ( - lim lz M (r) + \z,{r) + ^A 2 {r)Z,{r) ) , (8.8) 



which essentially differs from (17.91 ) by the absence of a finite contribution equiv- 
alent to \C[\]Zh{t) so that there is a problem in getting the correct Schott term. 
In fact, comparing with (T7.10I) . we see that the fundamental difference is that in 
Barut's approach the divering mass term is the 'Coulomb self-energy,' whereas in 
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our rigorous formalism it is a '5 2 -integral.' A similar situation arises in Dirac's 
approach [0, in which instead of calculating the self-interaction force the flow 
of energy-momentum out of a 3-surface surrounding the world-line is calculated, 
i.e., Eq. (11.81) . which implies 

QtiLir) = e 2 (-limlz M (r) + \a\t)Z»{t)), (8.9) 

in which there is no Z^t) contribution at all. 

Therefore, with self-interaction forces like (18.81) or (18.91) . there is no other 
possibility than to make one or more additional suppositions in order to get the 
Lorentz-Dirac equation of motion. In a way or another this amounts to modifying 
or supplementing Maxwell's theory and/or the notion of causality which leads to 
the retarded potentials, what is not necessary with the self-interaction force (17.101) 
in which the 5 2 -integral supplies the missing Z\i(r) contribution. 
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